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Abstract 

Self- avoiding walks and self-avoiding trails, two models of a polymer coil in dilute solution, have 
been shown to be governed by the same universality class. On the other hand, self-avoiding walks 
interacting via nearest-neighbour contacts (ISAW) and self-avoiding trails interacting via multiply- 
visited sites (ISAT) are two models of the coil-globule, or collapse transition of a polymer in dilute 
solution. On the square lattice it has been established numerically that the collapse transition of 
each model lies in a different universality class. 

The models differ in two substantial ways. They differ in the types of subsets of random walk 
configurations utilised (site self-avoidance versus bond self-avoidance) and in the type of attractive 
interaction. It is therefore of some interest to consider self-avoiding trails interacting via nearest 
neighbour attraction (INNSAT) in order to ascertain the source for the difference in the collapse 
universality class. Using the flatPERM algorithm, we have performed computer simulations of 
this model. We present numerical evidence that the singularity in the free energy of INNSAT at 
the collapse transition has a similar exponent to that of the ISAW model rather than the ISAT 
model. This would indicate that the type of interaction used in ISAW and ISAT is the source of 
the difference in universality class. 



1 Introduction 



The collapse transition of a polymer in a dilute solution has been a focus of study in lattice statistical 
mechanics for decades [TJ [2] . Any lattice model of a collapsing polymer has two key ingredients: an 
excluded volume effect expressing the impenetrability of monomers, and a short-range attractive force, 
which mimics the complex monomer-solvent interaction. When the effects of the excluded volume and 
the short-range attraction balance each other, the polymer undergoes a collapse transition which 
separates two distinct phases: a swollen and a collapsed phase. 

The canonical lattice model of the configurations of a polymer in solution has been the model 
of self-avoiding walks (SAW) where a random walk on a lattice is not allowed to visit a lattice site 
more than once. SAW display the desired excluded volume effect and are swollen in size relative to 
unrestricted random walks at the same length. A common way to introduce a short-range interaction 
is to assign a negative energy to each non-consecutive pair of monomers lying on neighbouring lattice 
sites, modelling an effective attractive force. This is the interacting self-avoiding walk (ISAW) model, 
which is the standard lattice model of polymer collapse using self-avoiding walks. 

The properties of lattice polymers are also related to those of magnetic systems near their critical 
point [3j. More precisely, lattice polymers are related to magnets with 0{n) symmetry in the formal 
limit of zero components (n — ¥ 0). This relation is of great importance, since it allows the application 
of the methods of statistical field theory to the study of polymer models, and the collapse transition 
can be understood as the tri-critical point of such systems [U IU [5] . 

The study of the critical properties of lattice polymers, and thus of O(n) models when we let n — > 0, 
in two dimensions has been ongoing over decades theoretically and numerically. Nienhuis in 1982 [6] 
was able to compute the critical exponents of free SAW by considering a model of non-intersecting 
loops on the hexagonal lattice, and in 1987 Duplantier and Saleur [7] were able to model the bond 
interaction introducing vacancies on the same lattice, obtaining a full set of critical exponents for 
the polymer collapse transition in the ISAW model. Their conjectured values for the exponents have 
been subsequently confirmed numerically by Prellberg and Owczarek in 1994 J5]. We will refer to (the 
universality class of) this critical point as the '0-point'. 

Introducing attractive interactions between bonds opens the doors to more complex possibilities. 
In the quest for a solvable 0(n) model on the square lattice, Blote and Nienhuis in 1989 |9j considered 
a lattice model, which includes energies for site-collisions and straight segments. For this model five 
critical branches are exactly known fTUj. [TTj 112] . One of these branches (named 'branch 0' in [§]) 
is similar to the #-point previously found by Duplantier and Saleur [7] (for this branch Batchelor 
[T3] obtained the exponents v = 4/7, 7 = 6/7), and two branches correspond to dense and dilute 
(SAW) polymers. The two remaining branches are, respectively, associated with a combination of 
Ising-like and 0{n) critical behaviour and with a new tri-critical point. This new tri-critical point 
(with exponents v = 12/23, 7 = 53/46) is another candidate for describing a collapsing polymer. 

A different model of a collapsing model can be constructed starting from self-avoiding trails. A 
self-avoiding trail (SAT) is a lattice walk configuration where the excluded volume is obtained by 
preventing the walk from visiting the same bond, rather than the same site, more than once. This is 
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a slightly weaker restriction, and SAW configurations are a proper subset of SAT configurations. The 
interacting version of self- avoiding trails (IS AT), customarily obtained by giving an energy to multiple 
visited sites, also presents a collapse transition. 

It is known that SAW and SAT share the same statistics in their high-temperature phase, but 
theoretical prediction and numerical evidence |14} [15] strongly suggests that the collapse transition 
of the ISAT model is in a different universality class to that of ISAW, although there is no clear 
understanding of why this would be the case. An interesting possible explanation of the ISAT collapse 
has been provided by Foster [16] in terms of the tri-critical point with exponents v = 12/23, 7 = 53/46 
derived from the 0{n) model. 

Although they both aim to describe the same physical system, the ISAW and ISAT models differ 
in both their geometrical properties and their interaction. Their different critical behaviour could be 
due to either one of these two differences. 

To investigate this further, we considered a mixed model where trails interact in the same way 
as SAW, that is, by a nearest-neighbour interaction. We shall call this new hybrid model Interacting 
Nearest-Neighbour Self- Avoiding Trails, or INNSAT. In the next section we review the ISAW and ISAT 
models describing the different behaviours that they have been shown to demonstrate. In Section [3] 
we formally introduce our new model INNSAT and the quantities of interest. We then describe our 
results in Section [4] and summarise our conclusions in the final section. 



2 ISAW and ISAT 

2.1 Interacting Self- Avoiding Walks (ISAW) 

Let us recall briefly the definition and main properties of the ISAW model. Consider the ensemble S n 
of self-avoiding walks (SAW) of length n, that is, of all lattice paths of n steps that can be formed on 
the square lattice such that they never visit the same site more than once. Given a SAW (p n 6 S n , we 
define a contact whenever there is a pair of sites that are neighbours on the lattice but not consecutive 
on the walk. We associate an energy — e c with each contact. Denoting by m c (ip n ) the number of 
contacts in cp n , the probability of <p n is given by 

e /3e c m c (ip n ) 



Z n (T) ' 

and the partition function Z n (T) is defined in the usual way as 



(2.1) 



Z n (T) = e^ m ^\ (2.2) 

where (5 is the inverse temperature l/fegT (kg is Boltzmann's constant). We define a Boltzmann 
weight (fugacity) uj c = exp(/3e c ). The finite- length reduced free energy is 

Kn(T) = — log Z n (T) (2.3) 
n 

and the thermodynamic limit is obtained by taking the limit of large n, i.e., 

k(T) = lim K n (T). (2.4) 
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As mentioned above, it is expected that there is a collapse phase transition at a temperature Tg, which 
is known as the #-point, characterised by a non-analyticity in k(T). 

The temperature Tg also separates regions of different finite-length scaling behaviour for fixed 
temperatures. Considering this finite-length scaling, for high temperatures (T > Tg) the excluded 
volume interaction is the dominant effect, and the behaviour is universally the same as for the non- 
interacting SAW problem: for large n, the mean squared end-to-end distance (or equivalently the 
radius of gyration) R\ and partition function Z n are expected to scale as 

R 2 n ~An 2v with v> 1/2 and (2.5) 
Z n ~Dv n n<-\ (2.6) 

respectively, where log// = k(T) and the exponents v and 7 are expected to be universal. In two 
dimensions it is well established [6] that v = 3/4 and 7 = 43/32 for T > Tg. The constants A and D 
are temperature dependent. 

Fixing the temperature precisely at the 0-point, T = Tg, Duplantier and Saleur [7] found v = 4/7 
and 7 = 8/7 in two dimensions. 

For low temperatures (T < Tg) it is accepted that the partition function is dominated by configura- 
tions that are internally dense, though not necessarily fully dense. The partition function should then 
scale differently from that at high temperatures, since a collapsed polymer should have a well-defined 
surface (and associated surface free energy) [17]. One expects in d dimensions large- n asymptotics of 
the form 

R 2 n ~An 2/d and (2.7) 
Z n r.D^f- 1),d n^-\ (2.8) 

with fj, s < 1. The constants A and D are temperature dependent. It is expected that the internal 
density smoothly goes to zero as the temperature is raised to the #-point. 

To explore the singularity in the free energy at the collapse point further, it is useful to consider 
the (reduced) internal energy and the specific heat, which are defined as 

u n (T) = and c n (T) = , (2.9) 

n n 

with limits 

U(T) = lim u n {T) and C(T) = lim c n (T). (2.10) 

n— >oo n— >oo 

When T — > Tg, the singular part of the specific heat behaves as 

C(T) ~ B \T 6 - Tp a , (2.11) 

where a < 1 for a second-order phase transition. If the transition is second-order, the singular part of 
the thermodynamic limit internal energy behaves as 

U (T) ~ B \Tg — T| 1_a (2.12) 
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when T — > Tq, and there is a jump in the internal energy at Tq if the transition is first-order (an 
effective value of a = 1). 

Tri-critical scaling |18] predicts that around the critical temperature, the finite-length scaling of 
the singular part of the specific- heat c n (T) obeys the following crossover scaling form 

c n (T)~n a( t> C{(T-T )n<t>), (2.13) 

when T — >• Tq and n — > oo, and that the exponents a and (j) are related via 

2-a=-. (2.14) 

9 

If one considers the peak of the finite-length specific heat it will behave as 

c p n eak ~n a(t> C(x max ), (2.15) 

where x max is the location of the maximum of the function C(x). 

The work of Duplantier and Saleur (1987) predicts the exponents for the #-point collapse are 

= ^ = 3/7 (2.16) 

and 

a = ag = -1/3. (2.17) 

It is important to observe that this implies that the specific heat does not diverge at the transition 
since the exponent 

a e( f) g = -1/7 « -0.14 (2.18) 

is negative. However, the peak values of the third derivative of the free energy with respect to 
temperature will diverge with positive exponent 

(1 + a e )<t>e = 2/7 « 0.28. (2.19) 
2.2 Interacting Self- Avoiding Trails (ISAT) 

The model of interacting trails on the square lattice is defined as follows. Consider the ensemble T n 
of self-avoiding trails (SAT) of length n, that is, of all lattice paths of n steps that can be formed on 
the square lattice such that they never visit the same bond more than once. Given a SAT ip n £ Tn, 
we associate an energy — et with each doubly visited site. Denoting by mt(i/j n ) the number of doubly 
visited sites in ip n , the probability of ip n is given by 

e /3etm t (4'„) 



zl SAT (TY 



(2.20) 



where we define the Boltzmann weight ujt = exp(/3et) and the partition function of the ISAT model is 
given by 

Z I n SAT {T)= Y, ^h Mn) - (2-21) 
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Previous work [141 [T5] on the square lattice has shown that there is a collapse transition at a tem- 
perature T = Tt with a strongly divergent specific heat, and the exponents have been estimated 
as 

(j) lT = 0.84(3) and a IT = 0.81(3) , (2.22) 
arising from a scaling of the peak value of the specific heat diverging with exponent 

IT a IT = 0.68(5) . (2.23) 

This result is a clear difference to the ISAW #-point described above where the singularity in the 
specific heat is convergent. Additionally, at T = Tt the finite-length scaling of the end-to-end distance 
was found to be consistent [H] with the form 

R 2 n {T) ~ An (Inn) 2 (2.24) 

as n — > oo. Again, this is quite different to the exponent v = 4/7 for the ISAW. 

Another important difference is that it has been recently observed [191 !i_20j that the low temperature 
phase is maximally dense on the triangular and square lattices. On the square lattice this implies that 
if one considers the proportion of the sites on the trail that are at lattice sites which are not doubly 
occupied via 

n — 2(mt) 

Pn = — , 2.25 

n 

then it is expected that 

p n — > as n — > oo. (2.26) 



3 The INNSAT model 

We define our new interacting model of trails (INNSAT) as follows. Consider the set of bond-avoiding 
paths T n as defined in the previous section. When two sites are adjacent on the lattice but not 
consecutive along the walk, so as not to be joined by any step of the walk, we again refer to this pair 
of sites as a nearest-neighbour contact and we give it a weight co c = e^ £c , analogously to the ISAW 
model. In Figure [I] a trail of length n = 21 with m c = 6 contacts is illustrated. 

Denoting by m c (ip n ) the number of contacts in ip n , the probability of ifj n is given by 



e Pe c m c (ip„) 



(3.1) 



where the partition function is 



Z^ T (T) = ^r c0M - (3-2) 

The intensive reduced internal energy and specific heat are as for ISAW: 

u n (T) = and c n (T) = . (3.3) 

n n 
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Figure 1: An example of INNSAT configuration with m c = 6, that is, there are 6 nearest-neighbour 
contacts illustrated via zigzag (red) lines. The trail can visit a site of the lattice twice by "touching" 
and by "crossing" itself. The number of doubly visited sites is rrit = 2. Note that there is no contact 
between the second and the seventh visited site of the walk, even though these are non-consecutive 
nearest-neighbour sites, as both sites are visited consecutively by a different segment of the trail. 

We shall also consider the proportion of the sites on the trail that are at lattice sites which are not 
doubly occupied via 

Pn= n ^ m ^, (3-4) 

n 

where wi; is the number of doubly visited sites as defined above. 

We shall see there is a collapse transition at a single temperature, T p , equivalently, fugacity oj v c . 
If the INNSAT model behaves in the same way as the ISAW model we would not expect the specific 
heat to diverge at this collapse point so we also need to introduce a quantity t n proportional to the 
third derivative of the free energy as 

fn = -3(m c )(m 2 c ) +2(m c ) 3 ^ 

n 

This quantity should have a singular part that behaves as 

t n (T) ~ n< Q+1 ^ T((T - T p )n^) , (3.6) 

around the collapse point, i.e., as T — > T p and n — > oo. Hence the peaks of this quantity t n will scale 
with exponent (a + 1)4>. 

4 Numerical Results 

We have first simulated the INNSAT model using the flatPERM algorithm [21] up to length n = 1000. 
With S n ~ 2.7 • 10 6 iterations, we collected 2.6 • 10 10 samples at the maximum length. Following [21], 
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we also measured the number of samples adjusted by the number of their independent growth steps, 
obtaining Sn — 1.6 • 10 s "effective samples". 

FlatPERM outputs an estimate W n mc of the total weight of the walks of length n at fixed values 
of m c . From the total weight one can access physical quantities over a broad range of temperatures 
through a simple weighted average, e.g. 



UJ 



w, 



n,m c 



(4.1) 



UJ',, 



We have begun by analysing the scaling of the specific heat by calculating the location of its peak 
= argmax^, c n (uj) and thereby evaluating (? n = c n (ujn)- We have also found the two peaks of t n (T) 
as a function of temperature: the peak values we denote tn^. 

It is clear that the specific heat peak is growing rather weakly as length increases. As the specific 
heat might converge, we obtain finite size estimates of the specific heat exponent by considering 



log 2 



C n/2 



C n/A 



(4.2) 



which should converge to a<p as n — > oo. Assuming corrections to scaling of n 3 / 7 as for ISAW, we 
find 

a NT <p NT = -0.15(5), (4.3) 

as shown on the left-hand side of Figure [2] This is consistent with a value of —1/7 ~ —0.14 for the 
#-point universality class. 
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Figure 2: Left: A plot of (4.2) designed to estimate the exponent a<p from the specific heat. Right: A 



plot of (4.4) designed to estimate the exponent (1 + a)4> from the derivative of the specific heat. The 



upper (green) points are values from the negative peak of t n while the lower (blue) points are from 
the positive peak of t n . 
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Furthermore when we analyse t n , which is proportional to the third derivative of the free energy, 
we find further consistent values of exponents. This derivative has two peaks tn^, and we obtain finite 
size estimates of the specific heat exponent from each of the peaks by considering 



log 2 



n/2j 



(4.4) 



which should converge to (1 + a)<f> as n 
estimate 



oo. From the plot on the right-hand side of Figure [2] we 



(a NT + 1)0 NT = 0.26(5). 



(4.5) 



There is still considerable curvature in the right-hand plot which would indicate a slightly larger value 
when higher order corrections to scaling are taken in account. Once again we conclude that the 
value of 2/7 tss 0.28 is consistent with our data. 
We therefore conjecture that 



and a NT = ag 7^ a r 



(4.6) 




Figure 3: A plot of the locations of the two peaks of t n (T) and of c n (T) against n~ 3 / 7 . Extrapolation 
for n — > 00 provides an estimate for the thermodynamic transition temperature cv p . The upper (green) 
curve is the location of negative tn~ peak, the lower (red) curve is the location of positive peak 
while the middle (blue) curve is the location of the peak of the specific heat. 

Using the locations of the peaks of the specific heat and its derivative and assuming finite-size 
correction of the order of n -3 / 7 , we have estimated the location of the collapse point as indicated in 
Figure [3j The magnitude of the difference between the three peak positions at finite length indicates 
the large corrections to scaling present in this problem at length 1000. We estimate 

uj p = 1.78(1). (4.7) 
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We next analysed the scaling of the end-to-end distance of the polymer. To obtain an estimate for 
the exponent v we considered a finite-size effective exponent 



,est 



(u) = I log 



(R 2 )n 



/2 



(4.8) 
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Figure 4: Estimates of the exponent v for INNSAT from finite-size estimates (4.8) as a function of 
temperature (left) and a zoom (right) on the crossing region. There are crossing points in the region 
of u between 1.78 and 1.80 giving corresponding estimates of v between 0.57 and 0.55. 

Using the data obtained from our flatPERM simulations, we plotted this quantity against temper- 
ature for various values of the length n, as shown in Figure [4j We find that the graphs for different 
values of n intersect around a particular value of the temperature. The location of this intersection 
point is a good estimator of the infinite-length critical temperature oj p and of the exponent v at the 
transition. We find a region of crossing points in uj between 1.78 and 1.80 with corresponding estimates 
of v between 0.57 and 0.55. While the obtained exponent estimate range is a little below the 6 point 
value of 4/7 ~ 0.571, we expect that our estimate is affected by strong corrections to scaling. For 
comparison, we simulated the ISAW model and performed an identical analysis, shown in Figure [5] 
The ISAW estimate of v ~ 0.56 at the collapse transition, while being closer to 4/7, also indicates the 
presence of strong corrections to scaling. 

Hence, we ran a thermal simulation at the estimated critical fugacity of 1.78 up to n = 8192. With 
5 = 5- 10 6 iterations, we obtained S n — 1.8 • 10 7 samples at the maximum length (corresponding to 
Sn ff ~ 1.2 • 10 5 "effective samples"). Figure 6 shows a log-log plot of the radius with respect to the 
size of the walk along with a linear best fit with slope 0.575. Of course, using u) = 1.77 orw = 1.79 
gives estimates of v that differ by 0.01 so the sensitivity to the location of the critical point dominates 
the error of our estimate. Correspondingly, if we assume that the INNSAT is in the ISAW universality 
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Figure 5: Estimates of the exponent v for ISAW from finite-size estimates (4.8) as a function of 
temperature (left) and a zoom (right) on the crossing region. The crossing point is at cj w 1.96 and 
v ~ 0.56. While the estimate for v at the crossing point is closer to the value of 4/7 than the one 
obtained for INNSAT, the difference provides evidence for the presence of fairly strong corrections-to- 
scaling. 

class we can use the value of v = 4/7 to better estimate the critical fugacity in the INNSAT model as 

uj p = 1.783(5). (4.9) 



4.1 Low-temperature phase 

There is strong evidence that the low-temperature phase of the ISAW model is a globular phase that 
is not fully dense, while for interacting trails the low-temperature phase is maximally dense. As 
discussed above, the trail fills the lattice asymptotically in a maximally dense phase, and the portion 
of steps not involved with doubly- visited sites should tend to zero as n — > oo. Following the analysis 
in [19] we measured the proportion p n of steps visiting the same site twice, and we plotted p n against 
n -1 / 2 at two different temperatures respectively above and below the critical temperature, as shown in 



Figure 4.1). We find that in both low-temperature and high-temperature phases only a small portion 
of the visited sites is visited again. The quantity p n stays close to 1 and, in particular, does not tend 
to zero in either region. This phenomenon can be understood by noting that the nearest-neighbour 
interaction makes energetically favourable for the walk to bounce away from an already visited site 
rather than to visit it again. This result implies that the low temperature phase is not fully dense, just 
as for the ISAW low temperature phase. This provides a consistent picture that the collapse transition 
in the INNSAT model is between high and low temperature phases similar to the ISAW, making the 
conclusion that the thermodynamic transition between them is similar to the 0-point a natural one. 
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Figure 6: Estimate of the exponent v from from finite-size estimates (4.8) obtained from a thermal 
simulation at oj = 1.78 up to length n = 8192. This estimate is sensitive to the estimate of the critical 
fugacity. A lower critical estimate would lead to a higher exponent estimate. 

5 Conclusions 

We have considered a new lattice model of polymer collapse in two dimensions which uses the inter- 
action type, nearest neighbour contacts, of the canonical lattice model (ISAW) but the configuration 
space of an alternative model (ISAT), being bond-avoiding walks, also known as self-avoiding trails. 
We find that the critical behaviour of the free energy, and its derivatives of our model (INNSAT) seems 
to align with the 0-point universality class of the ISAW model rather than that of the ISAT model, 
which has rather a strongly divergent specific heat not seen in either the ISAW or INNSAT models. 
The end-to-end distance scaling is broadly consistent with this conclusion, though estimates tend to 
be a little low. However, using the assumption that the INNSAT and ISAW models lie in the same 
universality class leads to a precise estimate of the critical fugacity as 1.783(5). Given the different 
universality classes for ISAT and INNSAT it may also be of interest to consider a model based upon 
self-avoiding trail configurations with both types of interaction. 
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